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ABSTRACT 

The complex structure of the hght curves of Swift GRBs has made the identification 
of breaks, and the interpretation of the blast wave caused by the burst, more difficult 
than in the pre-Swift era. We aim to identify breaks, which are possibly hidden, and 
to constrain the blast wave parameters; electron energy distribution, p, density profile 
of the circumburst medium, k, and the continued energy injection index, q. We do 
so by comparing the observed multi- wavelength light curves and X-ray spectra of our 
sample to the predictions of the blast wave model. We can successfully interpret all 
of the bursts in our sample of 10, except two, within this framework and we can 
estimate, with confidence, the electron energy distribution index for 6 of the sample. 
Furthermore we identify jet breaks in a number of the bursts. A statistical analysis of 
the distribution of p reveals that, even in the most conservative case of least scatter, 
the values are not consistent with a single, universal value. The values of k suggest 
that the circumburst density profiles are not drawn from only one of the constant 
density or wind-like media populations. 
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1 INTRODUCTION 

The afterglow emission of Gamma-Ray Bursts (GRBs) 
is generally well de scribe d by the blast wave model 



jRees fc MeszarosI [l993 : iMeszaros fc Reed 1 19971 : 
iMeszaros et al.l 1 19981 ). This model details the tempo- 
ral and spectral behaviour of the emission that is created 
by external shocks when a coUimated ultra-relativistic 
jet ploughs into the circumburst medium, driving a blast 
wave ahead of it. The level of coUimation, or jet opening 
angle, has important implications for the energetics of the 
underlying physical process, progenitor models, and the 
possible use of GRBs as standard candles. The signature 
of the coUimation, according to simple analytical models, 
is an achromatic temporal steepening or 'jet break' at 
approximately one day in an otherwise decaying, power-law 
light curve; from the t ime of this b reak, the jet opening 

angle can be estimated l|Rhoadslll997l ). 

Since the launch of the Swift satellite (|Gehrels et al.l 

I2OO4I ') it has become clear that this model for GRBs cannot, 
in its current form, explain the full complexity of observed 
light curve features and the lack of observed achromatic 
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temporal breaks. The unexpected f eatures detected, such 
as steep decays, plateau phases (e.g.. [ Tagliaferri et al. I I2OO5I : 
iNousek et"aLll2006l : [O'Brien et al.ll2006l') and a large number 
of X-ray flares (e.g.. iBurrows et al' 1 120071 : [Ch incarini et al.l 
I2OO7I : iFalcone et al.l l2007t ) have revealed the complexity of 
these sources up to about one day since the initial event, 
which is yet to be fully understood. These superimposed 
features also make it difficult to measure the underlying 
power-law features on which the blast wave model is based, 
and may lead to misinterpretations of the afterglows. 

Achromatic temporal breaks are not observed in the 
majority of bursts, up to weeks in a few bu rsts (e.g., 
iPanaitescu et al.ll2006l : H arrows fc Racusin|[2007l) . In some 
afterglows, a break is unobserved in both the X-ray and 
optical light curves, while in other bursts a break is ob- 
served in one regime but not the other (e.g., iLiang et al.l 
(2003). Previous compreh ensive studies of Swift lig ht curves 
in multiple wavelengths (|Liang et al.l |2007| . 120081 ) have at- 
tempted to identify achromatic temporal breaks due either 
to jet breaks, or the cessation of continued energy injection 
(i.e. end of the plateau phase), mainly as they concern the 
energy requirement of bursts. 

Identification of jet breaks is also important as a test 
of the blast wave physics as it represents the time at which 
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sideways spreading of the jet becomes important and the 
edges of the jet become visible. The post jet break phase 
also gives a direct measurement of the electron energy distri- 
bution, p, as the temporal decay, at all wavelengths, during 
this phase should have a value equal to that of p. How- 
ever, due to differences bet ween broken power- l aw fit s and 
the underlying parameters l|j6hannesson et al.ll2006l ). this 
method may systematically underestimate the true value of 
p. The value of p may also be measured from the spectral 
information of the optical or X-ray bands; the accuracy of 
this method is then dependent on the quality of the spectra 
used. A study of BeppoSAX bursts used broadband opti- 
cal to X-ray Spectral Energy Distributions (SEDs) to accu- 
rately constrain the value of p via the spectral and temporal 
parameters o f both bands, and a comparison to the blast 
wave model (|Starling et al.l l2008h . The same authors also 
use this method to estimate the density profile of the cir- 
cumburst medium, k, as the temporal slope may depend 
on this value as defined by p oc r~''; a constant density 
medium, fc = causing a more shallow decay than a wind 
like medium, k — 2. Similar interpretations of the tempo- 
ral and spectral information may also be used to estimate 
the effect, if any, that continued energy injection has on the 
light curves. Continued energy injection is used to explain 
the shallow decay phase observed in many Sw ift GRBs and 
is usu ally assumed to take the form oi E act'' (jNousek et al.l 
bOOd ). The additional energy can be due to slower shells of 
matter catching up with to blast wave, or a continued ac- 
tivity of the central engine. In either case, the addition of 
the extra energy to the blast wave causes the flux to drop 
off less rapidly than it would otherwise, and is observed as a 
shallower than expected decay. The cessation of the energy 
injection is marked by a temporal break to the normal decay 
phase as predicted by the blast wave model which implies, 
if it is due to continued activity of the central engine, that 
the engine must be active for up to about a day. Such late 
activity of the ce ntral engine has bee n observed in the form 
of episodic flares l|Curran et al.ll2008l l but such late, smooth 
and continuous activity has yet to be confirmed. 

In this paper we interpret a sample of afterglows of 
Swift GRBs, to constrain the blast wave parameters: elec- 
tron energy distribution, p, density profile of the circum- 
burst medium, k, and the continued energy injection index, 
q. In ^we introduce our sample and the method of reduc- 
tion, while in fJ3]we present the results of our temporal and 
spectral analysis. In 52] we interpret the results in the con- 
text of the blast wave model and discuss their implications 
in the overall context of GRB observations. We summarise 
our findings in 33 Throughout, we use the convention that 
a power-law flux is given as F oc t~°'v~^ where a is the 
temporal decay index and /3 is the spectral index. All uncer- 
tainties are quoted at the Icr confidence level. 



2 OBSERVATIONS AND ANALYSES 
2.1 Sample selection 

The bursts in our sample were chosen from an inspec- 
tion of previous li t eratu re on samples of temporal breaks 
jLiang et al. I l2007l . l2008l ). Bursts were also chosen from a 
comparison of the pre-reduced Swift X-ray Telescope (XRT, 



Table 1. Galactic absor ption, A^h jK alberla e t al.ll2005h . Galac- 



tic extinction, E, 



(B-V) 



dSchleeel et al . 1 998,) . and spectroscopic 



redshift, z, for the GRBs in our sample. 



GRB 




E(B-V) 


z 




xlO^O cm-2 




050730 


2.99 


0.051 


3.967" 


050801 


6.33 


0.096 


1.56+' 


050802 


1.86 


0.021 


1.71^= 


050820A 


4.41 


0.044 


2.6147'' 


051109A 


16.1 


0.190 


2.346= 


060124 


8.98 


0.135 


2.297^ 


060206 


0.889 


0.012 


4.0489 


060729 


4.49 


0.054 


0.54'' 


061121 


3.99 


0.046 


1.314' 


061126 


10.2 


0.182 


1.1588-'' 



t photometric redshift 

'iGhen et a.1.1 ll2005l'). 'jde Pasguale et al.l bOOTll. iFvnbo et al.l 



II2OO 5I). 'iLedoux et al.l ll2005li. HCenko et al.! 1I2OO6I). 
touimby et al.l ll2005h. ' Fvnbo et al.l ll2006h. iThone et al.l 
1I2OO6I) . iBloom et al.l ll20od '). - iPerlev et al.i 1 I2OO8I) . 



iBurrows et aP l2005l ) h ght curves in the on-line repository 
ijEvans et all l2007h up to the end of February 2008, and 
from the literature of optical data. 

Our sample (Table [ij consists of bursts with X-ray and 
optical light curves with good enough time coverage to allow 
for the underlying single power-law, or broken power-law, to 
be determined. The bursts are also well sampled enough in 
the X-ray and optical to constrain the spectral indices, /3x 
and /3opt. Bursts with published optical data were preferred 
over bursts with only GCN (Gamma-ray bursts Coordinates 
Network Circular) data which, by their very nature, are pre- 
liminary and hence difficult to compare with other data sets. 
However, in the cases of bursts with significant amounts 
of data from a single source, which we assume to be self- 
consistent, these data were included even if published only 
in GCNs. We did not confine our sample to bursts with clear 
breaks in either the X-ray or optical bands as we wanted to 
include the possibility of h idden or not very ob vious breaks, 
particularly in the X-ray l|Curran et al.l [20081 ') . or late, un- 
observed breaks. 



2.2 X-ray 

The XRT event data for our sample bursts were initially 
processed with the FTOOL, xrtpipeline (vO.11.4). Bad 
columns and pile-up were tested for and corrected for where 
necessary. Source and background spectra from the Photon 
Counting and W indowed Timing mode data (PC and WT; 
iHih et aL|[2003 ) were extracted using suitable extraction re- 
gions. The PC mode source extraction regions were circular 
(30 pixel radius), circular with a central exclusion radius 
for sequences suffering pile-up, or circular with exclusion 
radii to exclude nearby sources which would have contami- 
nated the spectra. WT mode source extraction regions were 
rectangular, and WT mode data did not suffer from pile- 
up during the time ranges for which we extracted spectra. 
Spectra were binned to have >20 photons per bin and the 
vOlO response matrices were used. 

Our light curve analyses are carried out on the pre- 
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Figure 1. The power-law fits to the XRT (upper data) and optical (lower data) light curves of each of the bursts in our sample. Data 
in grey were not used in the power-law fits. Optical hosts, where available, were included in the models. 
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Figure 1. Continued. 



reduced, XRT lig ht curves from the on-hne repository 
11 Evans et al These liRht curves were fit with single or 

broken power-laws (jBeuermann et aDll999l ) where the early, 
flaring data were eliminated from the fit. For bursts where 
there was a possible light curve break (Figure[T]), spectra 
were extracted pre-break (early) and post-break (late). The 
log mid-time of these spectra, foariy and tiato, are given in 
TableO If there was no possible break, early and late spectra 
were extracted for times before and after the approximate 
log mid-time of the afterglow, to test for any possible spec- 
tral change. The early and late spectra (0.3-10.0 keV) were 
fit both individually and simultaneously (early-|-late) with 
absorbed power-laws in Xspec 11.3.2 (Table [5]). statis- 
tics were used and the Galact ic values of Co l umn D ensity, 
A'h , were fixed to the values o f lKalberla et al. I l|2005l ) (Table 



2.3 Optical 

Optical photometric points in various bands were taken from 
the literature (Table [5)) and combined via a simultaneous 
temporal fit as detailed in lCurran" et all (1200781 ). This fitting 
allowed us to find the common temporal slope of the optical 



data and the colour differences between bands. Using these 
colours, the optical data were then shifted to a common mag- 
nitude and converted into an arbitrary, scaled flux to pro- 
duce joint optical light curves (Figure[T]). These light curves 
were fit with single or broken power-laws (Beucrmann et al.l 
1999), including host galaxy contributions where known (Ta- 
ble|4]). Early times at which the underlying behaviour was 
ambiguous, or flaring, were excluded from the fit. Where we 
deemed the errors to be unrealistically small and where they 
gave unacceptable values of x'^ >> 1, they were increased so 
as to give ~ 1, specifically in GRBs 050730 (x5), 050801 
(x2), and 050820 A (x2). 



2.4 Broadband SEDs 

The fit of optical colour differences and X-ray spectra de- 
scribed above were used to produce broadband spectral en- 
ergy distributions (SEDs) at early and late times for each 
burst in our sample. Optical observations were interpolated 
to the log mid-time of the X-ray spectra (toariy, iiatc), i.e., 
only bands with observations closest to that time were in- 
cluded in the SED. The optical SED points are co rrected 
for Galactic extinction, _E(s_v) (jSchlegel et al.|[l998l ). if not 



Testing the blast wave model with Swift GRBs 5 



Table 2. Optical references for the GRBs in our sample. 



GRB 


Filters 


Reference 


050730 


B,V,R, r',I,i',J,K 


1 


050801 


UVM2,UVWl,U,V,Rc 


2 


050802 


UVW2, UVM2, UVWl, U, B, V 


3 


050820A 


U,B,V,g,R,Ic,z' 


4 


051109A 


rR,CR 


5 


060124 


B,V,Rc,Ic 


6 


060206 


V,R,I,J,H,Ks 


7 


060729 


UVW2, UVM2, UVWl, U, B, V 


8 


061121 


UVW2, UVM2, UVWl, U, B, V, R 


9 


061126 


B,Rc,r',i' 


10 



llPandey et a l. (2006'!. 2 'Rvkoff ct al." l'2006''l:'dc Pasguale et al, 
20071'). S lOates ct al. ( 2007). 4 .Ccnko ct al. (2 06J, SlYost et al. 



2007| ). 6 ICurran et al.l 1 20061): iMisra et al.l (120071): Kann et 
in preparation ), 7_ Curran et al.l ll2007bl) : lAIatalo et al 



( 2006|) : IStanek et al.l( 2007): Wozni ak et al.l (120061). SiGrupe et al. 
20071). 9 iPage et al l"jl2'o07i) : jHalpern fc Armstrond (l2006allhl): 



iHalpern et al.l (l2006allbl) : iMelandri et al.l (|2006|) : lYost et al.l 
(|2006|) . lO lGomboc et aU (|2008|) . 



already done s o, and for lin e blanketing associated with the 
Lyman forest (|Madaulll995l). We used the simultaneous fit- 
ting method detailed in lStarling et al.l (|2007l ), where the op- 
tical to X-ray SED is fit in count-space, having the advan- 
tage that no model for the X-ray data need be assumed 
a priori. We fit the SEDs using models consisting of an ab- 
sorbed single or broken power-law (SPL or BPL) with slopes 
tied to Pi — 132 — 0. 5, as expected for a cooling break in th e 
blast wave model (|Meszaros et all Il998l : ISari et all 1 19981 ). 
We also use this data to calculate the spectral index between 
the optical at the lowest available frequency and the X-ray 
at 1.73 keV, /3opt-x (Table[S} to compare with the predicted 
values of an SPL, /3opt-x = /3x, or a BPL, /?opt-x = /3x-0.5. 

We use this broadband information in conjunction with 
the previously published SED analyses to estimate whether 
or not there was a break between the two regimes (TablelS)). 
When the of the SPL and BPL fits to the broadband 
SED are approximately 1, or the optical spectral index is 
consistent with the X-ray spectral index, we are unable to 
say if the broken power-law is a significantly better fit, so 
consider both cases. 



3 RESULTS 

Our results of the temporal (Figure[T| and X-ray spectral 
fits are summarised in Tables |3] and O and the details for 
individual bursts are discussed below. We defer the inter- 
pretation of these results to section |4] 



3.1 Individual bursts 

3.1.1 GRB 050730 

There is a clear break in both the optical and X-ray light 
curves at approximately the same time (~ 1 x 10* s). Fits 
of the X-ray light curve before this break are complicated 
by the lack of data and the possibility that the early flaring 
may continue to the break, hence the pre-break X-ray decay 
of Qx.i ~ 0.91 ± 0.07 should be treated with caution and 



as an upper limit. Similar deviations from power-law decay 
are observed in the optical but a greater time range of data 
before the break implies that the underlying slope here is 
less likely to deviate from the measured value. After the 
break the lack of optical data is reflected in the flt errors of 
the slope, which may be steeper than measured if there is a 
signiflcant flux contribution from a host (for which there is 
no estimate). 

There is an apparent X-ray spectral change across the 
break but the sparsity of optical data after the break makes 
this difficult to confirm. From our broadband fits we cannot 
make a conclusive statement on the existence of a spectral 
break but a single power-law offers an reasonable fit in agree- 
ment wi th our consistent X-r ay and optical to X-ray spectral 
indices. iPandev et al.l (|2006| ) find no evidence for a spectral 
break between the X-ray and optical, with an optical spec- 
tral index equal to that of the X-ray, at least before the 
break. 



3.1.2 GRB 050801 

A simultaneous break is observed in both bands and the op- 
tical light curve is well constrained both before and after the 
break. There is no detectable late break in the optical, as 
the last data point is consistent with the extrapolation from 
early times. The X-ray light curve is less well constrained due 
to few data points. Fitting the decaying X-ray data, after its 
peak at 290 s, with a single power-law, as opposed to a bro- 
ken power-law including the apparent rising segment, gives a 
significantly shallower decay of qx,2 = 1.16 ±0.03, in agree- 
ment with the optical. There is no measurable change in the 
X-ray spectral index over the break. Our SED fits are consis- 
tent with a SPL between the spectral regimes, in agreement 
with our optical to X-ray spectral indices being consi stent 
with the X-ray spectral i ndex. Both lRvkoff et alT (|2006| ) and 
Ide Pasguale et all (|2007l ) also find an X-ray to optical spec- 
tral index consistent with that of the X-ray, implying that 
there is no spectral break between the two regimes. 



3.1.3 GRB 050802 



In agreement with the results found bv lOates et al.l (|2007h . 
the X-ray light curve is well described by a broken power-law 
of ax,i = 0.70±0.04 andax,2 = 1.58±0.04, while the optical 
is described by a single power-law of ciopt,i = 0.81 ± 0.04. 
If a broken power-law is fit to the optical, the fit returns 
Q!opt,i = Q!opt,2, however the slopes and break time of the 
X-ray light curve may be forced, to increase the from 
0.89 to 1.47. This may be an acceptable fit since we cannot 
apply a correction for a host galaxy as it was not observed. 

We find no measurable change in the X-ray spectral 
index over the temporal break. From our broadband SED 
fits it is not possible to rule out a spectral break, though 
neither is one required. This is in agreement with an X-ray 
to optical spectral in dex consisten t with the X-ray spectral 
index as both we and lOates et al.1 (|2007D find. 

3.1.4 GRB 050820 A 

After an initial optical rise and X-ray fiaring, both light 
curves are well fit by broken power-laws with consistent 
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Table 3. F-test probabilities, Fprobi for the fits of the early time, toarlyi a-iid late time, tiatci broadband SEDs and whether they are best 
described by a single or a broken power-law (SPL or BPL). 



GRB 


^carly 
S 




-^prob 


*latc 
S 




-^prob 




050730 


7.03 X 


10^ 


0.38 


8.77 X 


10* 


1.0 


SPL/BPLt 


050801 


161 




0.20 








SPL/- 


050802 


1.18 X 


10^ 


0.01 


8.88 X 


10* 


8 X 10-* 


SPL/BPL 


050820A 


1.66 X 


10"* 


3 X 10"" 


7.59 X 




0.63 


-/BPL 


051109A 


8.96 X 


10^ 


1.0 


2.80 X 


10^ 


1.0 


SPL/BPLt 


060124 


3.19 X 


10-* 


3 X 10-21 


2.61 X 


10^ 


9 X 10-9 


-/BPL 


060206 


1.00 X 


10* 


1.9 X 10~^ 


8.20 X 


10* 


0.63 


SPL/BPL 


060729 


2.09 X 


10* 


9 X 10"" 


1.03 X 


10** 


3 X 10-5 


-/BPL 


061121 


9.66 X 


10^ 


6 X 10"^ 


2.63 X 


10^ 


0.02 


-/BPL 


061126 


2.59 X 


10^ 


6 X 10-^ 


1.44 X 


10^ 


1.0 


-/BPL 



t The of the SPL and BPL are approximately 1 so wc cannot state if the improvement is significant, despite the F-tcst. 

Table 4. Temporal fits of the light curves in our sample; pre-break temporal index, break time in seconds and post-break index for the 
X-ray and the optical light curves. 



GRB 


ax,l 


^X, break 

s 


Q!X,2 






^opt, break 
S 


a 


opt, 2 


050730 


0.91 ± 0.07 


(1.06 ±0.03) X 10* 


2.56 ± 


0.04 


0.55 ± 0.03 


1.4tJJ xlO* 


1. 


■ '"-0.12 


050801 


-0.6 ± 0.4 


290t|0 


1.39 ± 


0.04 


0.13 ± 0.03 


230 ± 10 


1, 


.16 ± 0.01 


050802 


0.70 ± 0.04 


(7.8 ±0.1) X 10^ 


1.58 ± 


0.04 


0.81 ± 0.04 








050820A 


1.17 ± 0.02 


(5.9 ± 1.6) X 10^ 


1.72 ± 


0.15 


0.91 ± 0.01 


(5.5 ± 1.0) X 10^ 


1, 


.60 ± 0.12 


051109A 


1.07 ± 0.03 


(8 ± 4) X 10* 


1.36 ± 


0.05 


0.66 ± 0.01 


(5.7 ±0.7) X 10* 


1, 


.54 ± 0.13 


060124 


1.01 ± 0.08 


(5.8 ± 1.0) X 10* 


1.45 ± 


0.03 


0.83 ± 0.01 


l.Otll, X 10' 


1, 


.35 ± 0.02 


060206 


1.04 ± 0.10 


2.2lg:° X 10* 


1.40 ± 


0.07 


1.022 ± 0.006 


(6.23 ±0.05) X 10* 


1, 


.99 ± 0.05 


060729" 


0.23 ± 0.02 


(7.4 ±0.3) X 10* 


1.28 ± 


0.02 


-0.26 ± 0.02 


(4.2 ±0.1) X 10* 


1, 


,28 ± 0.02 


061121'' 


0.97 ± 0.02 


(2.9 ±0.4) X 10* 


1.55 ± 


0.03 


0.96 ± 0.02 


2.9 xlO* (fixed) 


1, 


.55 (fixed) 


061126 


1.29 ± 0.01 








0.91 ± 0.01 


(1.33 ±0.08) X 10^ 


1, 


.9 ± 0.1 



" Fitting a doubly broken power-law to the X-ray, as plotted, also reveals a late break at (1.4 ± 0.4) X 10'' s to ax, 3 = 1-70 ± 0.09. 
'' Fitting a doubly broken power-law to the X-ray, as plotted, also reveals an early break at 2000 ± 400 s from Qx,o = 0.32 ± 0.04. 



break times and post break slopes. According to lCenko et al.l 
l|2006ll a host contribution may affect the last optical (HST) 
points at ~ 3 X 10*' s so the fit decay and break time should 
possibly be treated as lower limits. We find no need to fit a 
decaying segment from 500 - 5000 s, as they have done, as 
an extrapolation of the later data back to this time does not 
support a break. 

Our broadband SED fits and our derived X-ray and 
optical to X-ray spectral indices support a brea k between the 
two sp ectral regimes. This is in agreement with lCenko et all 
l|2006l ) who find an optical spectral index, (3opt ~ 0.67±0.08, 
3.9a shallower than the X-ray spectral index of /3x = 1.00 ± 
0.03. 



3.1.5 GRB 05 11 09 A 

There is a clear break in the optical light curve and even 
though it is not obvious, an F-test probability of 5 x 10"^, 
supports a break at a consistent time in the X-rays, with 
a confidence of almost 5a. Fixing the X-ray break time to 
the optical break time only affects the decay indices within 
uncertainties. From our broadband SED, as well as our X- 
ray and optical to X-ray spectral indices, we find both single 
and broken power-laws offer acceptable fits. This is in agree- 



ment with lYost et al.l (|2007[ ) who find an X-ray to optical 
spectral index of ~ 0.7 ± 0.2, only ~ 1.5o- lower than the 
X-ray spectral index of /3x = 1.04 ± 0.04. 

3.1.6 GRB 060124 

After the initial X-ray fiaring, when the optical is not sam- 
pled well enough to be able to judge its behaviour, both 
afterglows are well described as broken power-laws with 
consistent break times and post-break slopes of Qopt.i ~ 
1.35±0.02 and ax,i = 1.45±0.03. We find no X-ray spectral 
change over this break. Our broadband SED conclusively re- 
quires a spectral break between the two regimes, consistent 
with the ambiguous difference between the X-ray and optical 
to X-ray spectral indices. 

3.1.7 GRB 060206 

After initial rising in the optical and flaring behaviour in 
the X-ray light curves, both decay as power-laws. The opti- 
cal afterglow is clearly a broken power-law, while the X-ray 
is well fit by a single power-law of a = 1.28 ± 0.02. How- 
ever, the X-ray can also be fit by a broken power-law if the 
data past 1 x 10® s are eliminated. This is valid as there 
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Table 5. X-ray spectral fits for the GRBs in our sample; spectral indices, /3x, and host extinctions, A^Hi for early and late time (as 
defined in Table|3]l X-ray spectra (see section [2]2]l . The difference between /3x, early and /3x,latoi A/3xi is given as the number of a, where 
a is the quadrature sum of the errors. The optical to X-ray spectral index, /3opt-Xi is also given. 



GRB 


/3x, early 


-'Vh, early 

xlO'^^ cm^^ 


/^X.latc 


-^'^H.lato 

xlO^^ cm~2 


A/3x 

cr 


/5x,early+late 


-'Vn.early-l-late 

X lO-^-^ cm~2 


/3opt-X 


050730 


0.56 ± 0.03 


1.1 ± 0.2 


0.73 ± 0.03 


0.90 ± 0.3 


4.0 


0.63 ± 0.02 


1.08 ± 0.17 


0.85±0.06 


050801 


0.9 ± 0.2 


< 1.0 


0.79 ± 0.10 


< 4.25 


0.5 


0.82 ± 0.10 


< 0.24 


1.08±0.16 


050802 


0.75 ± 0.04 


0.28 ± 0.06 


0.89 ± 0.04 


0.30 ± 0.07 


2.5 


0.81 ± 0.03 


0.29 ± 0.05 


0.61±0.10 


050820A 


1.03 ± 0.04 


0.80 ± 0.12 


0.90 ± 0.05 


0.29 ± 0.17 


2.0 


1.00 ± 0.03 


0.68 ± 0.10 


0.37±0.15 


051109A 


1.04 ± 0.04 


1.1 ± 0.2 


1.05 ± 0.09 


1.0 ± 0.45 


0.1 


1.04 ± 0.04 


1.1 ± 0.2 


0.73±0.13 


060124 


1.02 ± 0.04 


0.91 ± 0.12 


1.00 ± 0.04 


0.90 ± 0.13 


0.4 


1.01 ± 0.03 


0.91 ± 0.10 


0.74±0.10 


060206 


1.26 ± 0.08 


1.2 ± 0.4 


0.92 ± 0.13 


1.3 ± 0.7 


2.2 


1.17 ± 0.08 


1.2 ± 0.4 


0.90±0.08 


060729 


1.11 ± 0.03 


0.16 ± 0.01 


1.12 ± 0.02 


0.25 ± 0.01 


0.3 


1.11 ± 0.02 


0.20 ± 0.1 


0.92±0.08 


061121" 


1.01 ± 0.05 


0.79 ± 0.09 


0.83 ± 0.04 


0.68 ± 0.10 


2.8 


0.94 ± 0.04 


0.74 ± 0.08 


0.59±0.21 


061126 


0.84 ± 0.06 


0.68 ± 0.13 


0.95 ± 0.03 


0.69 ± 0.08 


1.6 


0.91 ± 0.04 


0.69 ± 0.08 


0.50±0.29 


' For the first decay segment, before 2000 s 


wc find /3x = 0.9 


9 ± 0.07 and A^h 


= (0.78 ±0.13) X 10^2 cm 


-2 





is signific ant contamination fr om a nearby source as dis- 
cu ssed bvlCur ran ct al. ('2007b") and shown in the reduction 
of iButler fc"Ko cevslci (,2007 ). The post- break X -ray tempo- 
ral index should hence be treated as a lower limit. Neither 
the difference between the X-ray and optical to X-ray spec- 
tral indices nor the SED are conclusive as to whether a break 
is required but given that the F-test probability corresponds 
to a le ss than 3ff certainty we favour here, as we have previ- 
ously l|Curran et al.ll2007bh . a single power-law. Neither do 
we find any change in X-ray spectral index over the break 
greater than 3cr. 

3.1.8 GRB 060729 

Both light curves follow parallel slopes after an almost si- 
multaneous break at ~ 6 x 10* s which shows no X-ray 
spectral change. Before this break, from 5000 s, the X-ray 
light curve is decaying with ax,i = 0.23 ± 0.02, though this 
is somewhat dependent on the time one chooses to start 
the fit. The optical flux over this time is increasing with 
Qopt,i = —0.26 ± 0.02, though if the earlier optical data 
are included, this is reduced to Qopt.i = — O.OS^Q Qg, leaving 
the post-break slope unchanged, within uncertainties, and 
causing the break time to increase but not significantly. Fit- 
ting a doubly broken power-law to the X-ray light curve, 
we also find a second break in the X-ray light curve at 
(1.4 ± 0.4) X 10" s to ax, 3 = 1-70 ± 0.09. An F-test compari- 
son of the of the doubly broken power-law and the broken 
power-law give a highly significant probability of 2.5 x 10~^, 
which is much greater than 5(t. The optical light curve does 
not sample this time so is unable to confirm or refute the 
break, which was n ot identified in the light curve analysis of 
iGru pe et all l|2007h . 

Our SED analysis finds a BPL is favoured over a SPL 
and hence there is likely a spectral break between the two 
regimes; in direct contrast to the parallel nature of the op- 
tical and X-ray light curves. However the early time SED, 
where the spectral break is most pronounced, is taken before 
the temporal break at ~ 6 x 10* s, where the light curves are 
not parallel; the late time SED, d uring the parallel p hase, 
supports a BPL with less certainty. iGrupe et al.l l|2007l ) find 
a constant optical to X-ray spectral index of ~ 0.8 ± 0.05, 
compared to our X-ray spectral index of /3x = 1-11 ± 0.02 



and optical to X-ray spectral index of /3opt-x = 0.92 ± 0.08, 
all of which are consistent with a SPL. Given this we cannot 
conclusively support the BPL nature if the SED. 

3.1.9 GRB 061121 

The X-ray light curve displays a clear steepening from early 
to late times which is modelled by a doubly broken power- 
law. This a significantly better fit than a broken power-law 
with a break at ~ 9 x 10^ s (F-test probability ~ 10"^°). 
After a rise to 212 s, the optical light curve appears to decay 
with a single power-law of aopt.i = 0.92 ± 0.021. However , 
when the effects of the host galaxy, which [Page et all (|2007l ) 
estimate to be ~ 2.5^Jy in the R band, are taken into ac- 
count, a break is detectable. There is not enough data after 
the break to constrain the optical break time or temporal 
decay so they are set at the X-ray values, which gives an 
acceptable fit. 

As the light curve displays 3 power-law segments, X- 
ray spectra were extracted for each (Table[5]). The first two 
segments are in excellent agreement with each other and 
in marginal agreement (2.8cr) with the last segment. The 
first segment, and hence the second also, may be affected by 
flux from the early prompt emission. We flnd a significant 
improvement in the description of the optical to X-ray SED, 
when a broken power-law is used in favour of a single power 
law, in agreement wit h the optical to X -ray spectral index 
as weU as the SED of iPage et all (|2007h . 

3.1.10 GRB061126 

Before the effect of the host is taken into account, the op- 
tical light curve is well described by a single power law of 
Qopt = 0.930 ± 0.005. When the host is modelled, a break 
from Oopt.i = 0.91 ±0.01 to aopt,2 = 1.9±0.1 becomes clear. 
The X-ray light curve does not seem to be consistent with 
a break since when a broken power-law, with a break time 
equal to that of the optical, is flt we find ax.i = ax, 2. How- 
ever, forcing the post-break slope to equal that of the optical 
only increases the from 1.35 to 1.60 and while there is no 
obvious source of contamination to the X-ray light curve, we 
cannot rule out a break. We find in our broadband SED fit, 
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as do lGomboc et al.l l|2008l ). that the optical to X-ray SED 
is best fit by a broken power-law at least before the break. 
The SED after the break and the optical to X-ray spectral 
index are consistent with this, though inconclusive. 



4 DISCUSSION 

We use the blast wave model (|Rees fc Meszarod 1 19921 : 
iMeszaros et al.|[l998l ) to describe the temporal and spectral 
properties of the GRB afterglow emission. We assume on- 
axis viewing, a standard jet structure and no evolution of the 
microphysical parameters. The relations between the tempo- 
ral and spectral indices and the blast wave p arameters that 
we us e are summarised in Tab le[6](see also e.g. lStarling et all 
l2008l and lNousek et al.ll2006l ). Our general method was to 
estimate the value of the electron energy distribution index, 
p, from the X-ray spectral index (section l4.1|l and use this to 
calculate the predicted values of temporal decay (which are 
sub-scripted with p, section IT^ . We derive p from the spec- 
tral index as opposed to the temporal index since for a given 
spectral index there are only two possible values of p, while 
for a given temporal index there are multiple possible values. 
Spectral slopes are dependent only on p and the position of 
the cooling break. Temporal indices, a, are dependent on p, 
the position of the cooling break, the circumburst density 
profile, k, and on possible continued energy injection (Ta- 
ble|S]). We assume that continued energy injection takes the 
form oi E (xf where q is dependent on the shallowing effect 
of the injection, Ag = ap — aobsorvcd, and on the value of p 
jNousek et aLll2006 !: though they define q as E oc t''+^). This 
injection period from ti to tf increases the energy of the af- 
terglow by a factor of / = {tf /ti)'' . Temporal indices are also 
pron e to being incorrectly est imated from broken power-law 
fits (|J6hannesson et al.ll20od) . underestimating post- break 
indices. 

4.1 Derivation of p 

The values of the electron energy distribution index, p, are 
derived from the X-ray spectral index of the afterglow, /3x 
(Table [SJ, following the relations in Table [S] As most of the 
sample show no evolution of X-ray spectral index, we used 
the average spectral index, /3cariy+iato in all cases except one. 
For GRB 050730 there was an apparent spectral evolution so 
we used only the late time spectral index, /3iato, since this 
was less likely to be contaminated by flux from the prompt 
emission; in this case the late spectrum better overlaps with 
the time range where we apply the blast wave model to the 
light curve so it is more proper to use it alone. 

For a given value of X-ray spectral index, there are two 
possible values of p (TableO depending on whether the cool- 
ing break, Uc, is below (p — 2/3) or above (p — 2f3 + 1) the 
X-ray frequencies, i/x. If the optical to X-ray SED does not 
display a break then the cooling break can either be above 
the X-ray regime or below the optical regime and the blast 
wave predictions of each p are compared to the observed 
temporal slopes to discern which is correct. If the SED re- 
quires a broken power-law it most likely implies that a cool- 
ing break lies between the two regimes and is below the X- 
ray regime. A cooling break requires, or must be consistent 
with, a difference between the spectral slopes of A/3 — 0.5. 



However, a break between the two regimes does not neces- 
sarily imply a cooling break; it may be due to the fact that 
each regime has a different spectral index since they are orig- 
inating from different emis sion regions, as is cl aimed to be 
the case for GRB080319B (|Racusin et al.ll2008l ). A spectral 
break due to two different emission regions does not have a 
predictable difference between slopes but for this interpreta- 
tion to work, one must be able to explain why each emission 
region is only visible in one regime. A cooling break is a more 
likely explanation in the majority of cases but a comparison 
of the blast wave predictions of each p with the observed 
light curves is required. 

4.1.1 Broadband versus X-ray 

We use the X-ray spectral index to derive the value of p 
because the value of /3x is relatively unaffected by environ- 
mental unknowns. The main unknown in the X-ray spectral 
fits is that of the absorption due to the column density in the 
host galaxy; this is easily fit and only affects the low energies 
of the Swift XRT spectra. X-ray spectra may also suffer con- 
tamination from nearby sources but this can be minimised 
through careful extraction of spectra. The optical spectral 
index is dependent on, and suffers somewhat of a degener- 
acy with, extinction in the host galaxy. The extinction af- 
fects the whole range of the optical wavelengths somewhat, 
while absorption only significantly affects less than about 
half a decade of the X-ray range. In fact, as the redshift in- 
creases, extinction has a greater effect on the optical bands, 
while absorption is shifted out of the X-ray range. The opti- 
cal data ideally spans from the near infrared K band to the 
ultraviolet U band, offering approximately a decade of spec- 
tral range with which to constrain the spectral index. This 
is a range similar to that offered by the XRT (0.3-10 keV) 
but is rarely reached and, given the narrow range of com- 
monly observed bands (UBV or VRI), may not be able to 
constrain the spectral index very well. 

Broadband, optical to X-ray, spectra are used: to elimi- 
nate some of the difficulties associated with the optical data; 
increase the fitting range; and reduce the number of free 
parameters by tying the optical and X-ray spectral indices 
to each other. This has proved successful in a number of 
cases and especially for BeppoSAX (~0.1-10keV) bursts 
(|Starling et al.l l2008l where the quality of the X-ray spectra 
did not always allow the spectral index to be sufficiently 
constrained without the optical data; t he mean spectral 
index error of BeppoSAX data is 0.30 l|de Pasguale et aU 
; 20081 . Swift XRT allows for more accurate X-ray spectral 
indices with the m edian error on XRT PC spectra being 0.11 
lEvans et al.ll200^ but this overestimates the error on many 
bursts of interest as Swift detects many more faint bursts. 
In our sample the average spectral index error is only 0.043 
or less than 5% of the average spectral index of 0.94, from 
our early-(-late spectra. 

While the spectral index should only be constrained fur- 
ther in a broadband fit, the difficulties associated with the 
combination of data from different telescopes, may intro- 
duce errors and incorrect estimates of p. There may also be 
an unknown or uncertain contribution from a host galaxy, or 
there may be undocumented errors on the zero-points used 
to convert from magnitude to fiux. It is also worth noting 
that data used for the optical SEDs, if taken from liter- 
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Table 6. A summary of the relationships for temporal and spectral indices, a and /3, for a given value of electron energy distribution 
inde x, p, in the diff e rent r egimes (from e.g. .Starling et al...2008- . assuming p > 2). The shallowing effect, Aa, of continued energy injection 
from lNousek et al.l 1I2OO6I') . 







a 




Aa 




/3 


Vc,m < I' (k) 




(3p - 2)/4 


q{2 


+ p)/4 




p/2 


I^c,in < I' (jet) 




P 








p/2 


I'm < 1^ < l^c (k) 




12{p-l)-fc(3p-5) 
4(4-fc) 


q{3- 


fc + p)/4 


(P 


-l)/2 


<v < Vc {k = 


0) 


3(p-l)/4 


q(3 


+ p)/4 


iP 


-l)/2 


f m < < I'c (k = 


2) 


(3p-l)/4 


9(1 


+ p)/4 


(P 


-l)/2 


< U < Uc (jet) 




P 






(P 


-l)/2 



Table 7. Possible values of p, derived from the X-ray spectral 
index of the afterglow, /3x ■ 



GRB 






C 


050730 


1.46 


± 0.06 


2.46 




0.06 


050801 


1.64 


± 0.20 


2.64 


± 


0.20 


050802 


1.62 


± 0.06 


2.62 


± 


0.06 


050820A 


2.00 


± 0.06 


3.00 


± 


0.06 


051109A 


2.08 


± 0.08 


3.08 


± 


0.08 


060124 


2.02 


± 0.06 


3.02 


± 


0.06 


060206 


2.34 


± 0.16 


3.34 


± 


0.16 


060729 


2.22 


± 0.04 


3.22 


± 


0.04 


061121 


1.88 


± 0.08 


2.88 


± 


0.08 


061126 


1.82 


± 0.08 


2.82 


± 


0.08 



ature, may have to be converted into different quantities 
for fitting and these conversions can add significant errors. 
While broadband SEDs may offer excellent constraints on 
optical spectral indices, extinction and absorption, they are 
not generally necessary to constrain the X-ray spectral in- 
dex and hence the possible values of p. They are however 
invaluable in the identification of spectral breaks between 
the two regimes which is necessary in the interpretation of 
a GRB within the blast wave model and the correct choice 
of p. 

4.2 Blast wave interpretation 

Here we compare the derived values of p and the blast 
wave model predictions to the observed temporal and spec- 
tral properties to find the most likely parameters (Table|8]). 
Given the complexity and ambiguity of some of the tempo- 
ral results summarised in Table |3J we advise reference to the 
text of section [3] also. 

4.2.1 GRB 050730 

We find that p — 2.46 ± 0.06 for i^x < Vc and unconstrained 
k and g, though this in not unequivocal. The steep X-ray 
temporal decay of ax, 2 = 2.56 ± 0.04 in this burst can only 
be consistent with a post jet break decay of ax,p = Qjct = p 
given p = 2.46 ± 0.06 derived from the X-ray spectral in- 
dex, assuming vx < Vc. The optical slope while shallower is 
also consistent, within 2a, with a post jet break decay equal 
to p. This implies that parallel, pre-break slopes of between 
1.10 and 1.60, depending on k, or shallower given energy 
injection are required for both regimes in disagreement with 



the 4(7 difference between the observed slopes. Since we can- 
not confirm a possible steepening contribution to the X-ray 
light curve from fiares, the interpretation of this burst and 
the derived value of p is ambiguous. 

4.2.2 GRB 050801 

We find that p = 2.64±0.20 for i^x < Vc and k = -0.7± 1.4, 
q — 0.65 ±0.15. If we adopt the value of post-break temporal 
X-ray slope from the single power-law, as opposed to that 
from the broken power-law, both regimes are parallel after 
the break with ax ~ aopt.i = 1.16±0.01, in agreement with 
the lack of a break in the SED. This implies vx < Vc and 
p = 2.64 ± 0.20, which predicts a slope of ax,p = aopt,p = 
1.23±0.09 for fc = 0, in agreement with the observed values. 
Calculating k from the temporal slope and value of p we get 
k = -0.7±1.4. The break of Aa = 1.03±0.03 observed in the 
optical, could then be due to energy injection at early times 
of q = 0.65 ± 0.15, increasing the energy of the afterglow 
over th observed period of injection by a factor of / ~ 5. 
The behaviour of the X-ray at that time is unclear, though 
it does not rule out a break. 

4.2.3 GRB 050802 

We find that p = 2.62 ± 0.06 for vx < Vc and fc = 1.7 ± 0.3, 
q = 0.90 ± 0.13. The post-break X-ray slope of ax, 2 = 
1.58 ± 0.04 is consistent p = 2.62 ± 0.06 {vx < i^c) and a 
circumburst density profile of k = 1.7±0.3. Since it is not 
certain whether there is a break in the SED, the position- 
ing of the X-ray frequency below the cooling frequency is 
valid. Given that we cannot rule out a temporal break in 
the optical light curve, it is also possible that aopt = ax as 
required. The X-ray break of Aa = 0.88 ± 0.06, consistent 
with the optical, would then be due to energy injection of 
q = 0.9 ± 0.13 which increases the energy of the afterglow 
over the observed period of injection by a factor of / ~ 10. 

If there is not a break in the optical, then a possible 
explanation may be that the two regimes are originating 
from different emission regions as proposed bv lOates et al.l 
(2007), though this raises the question as to why only one 
regime is observed from each emission region. 

4.2.4 GRB 05 0820 A 

We do not find a confident blast wave interpretation consis- 
tent with the light curves of this burst. The inconsistency of 
the optical and X-ray slopes, in agreement with the broken 
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power- law SED, suggests that the coohng frequency is be- 
tween the two regimes. This would imply that p = 2.00±0.06 
and ax,p = 1.00±0.03, 4.7a shallower than the observed pre- 
break value of ax,i = 1.17±0.02. The optical, which is below 
the cooling frequency, is explained by a circumburst medium 
of density profile, A: = 1.0 ± 0.3. The post-break slopes are 
consistent with the predicted ajet = p = 2.00 ± 0.06 of this 
interpretation at the Scr level and possibly better given that 
the optical may be affected by a unquantified host contribu- 
tion. Given the 4.70- difference between the predicted and 
observed X-ray slopes before the break we cannot confi- 
dently support this interpretation. 

4.2.5 GRB051109A 

We find that p = 2.08±0.08 for < px and k = -1.7±0.8, 
g = 0. In this burst it is unclear whether or not there is a 
cooling break between the two regimes so we have considered 
both cases. li vx > Vc then p = 2.08 ± 0.08 implies an X-ray 
temporal slope of ax,p = 1.06 ±0.04, in agreement with the 
X-ray slope before the break. The optical is then below the 
cooling break and p implies Oopt.p = 0.81 ± 0.03 for k = 0, 
which is greater than the observed value Oopt = 0.66 ± 0.01 
by 4.7(7; in fact the observed temporal slope requires k = 
— 1.7 ± 0.8. The break, which is consistent in optical and 
X-rays is then interpreted as a jet break, which should go 
as Ojet = p = 2.08 ± 0.08. The observed value of a ~ 1.4 
is shallower than expected but may be rolling over to the 
predicted value. 

Alternatively, if < Vc then p = 3.08 ± 0.08 implies 
temporal slopes of ax,p = Oopt.p = 1.56 ± 0.03 with = 0, 
comparing well with the observed post-break slopes of both. 
The break can then interpreted as the cessation of continued 
energy injection. The difference between the observed pre- 
break indices and the values predicted by the blast wave are, 
on average, Aa ~ 0.7, which corresponds to g ~ 0.6. How- 
ever, the 13(7 difference between these slopes, which should 
be parallel, is difficult to explain, making our first explanar- 
tion more likely. 

1.2.6 GRB 060124 

We find that p = 2.02 ± 0.06 for ly^ < vx and A; = 0.0 ± 0.4, 
g = 0. The broken power-law SED and unequal pre-break 
decays require that the cooling break lies between the optical 
and X-ray frequencies, which implies p = 2.02 ± 0.06. The 
predicted slopes of ax,p = 1.02 ± 0.03 and aopt.p = 0.77 ± 
0.02 (k = 0.0 ± 0.4) are in agreement with the observed 
values before the break. The observed post-break values of 
~ 1.4 arc shallower than the predicted value of Ojct = p, 
but may be rolling over to the asymptotic value, or may be 
shallower than expected due to a structured jet as suggested 
by Kann et al. (in prepaxation) . 

4.2.7 GRB 060206 

We find that either p = 2.34 ± 0.16 for I'c < vx with q = Q 
and k unconstrained, or, p = 3.34 ± 0.16 for i>x < i^c and 
k = 1.3±0.7, q = 0.8±0.2; the data are ambiguous. The sin- 
gle power-law SED from optical to X-rays and the fact that 
both have equal temporal decay indices before the break. 



imply that both regimes are either above, or below, the 
cooling frequency. In the first case p = 2.34 ± 0.16, imply- 
ing Op = 1.26 ± 0.09 before the break, which is slightly 
steeper than, but within 3a of the observed slopes which 
are qx,i = 1.04 ± 0.10 and aopt.i = 1.022 ± 0.006. The pre- 
dicted post-break slope for a jet break, ajet = p, could agree 
with the optical slope rolling over to an asymptotic value 
and with the 'hidden' X-ray break. 

In the second case, the observed bands are below the 
cooling frequency and p — 3.34 ± 0.16, implying Op is be- 
tween 1.76 ± 0.08 (fc = 0) and 2.26 ± 0.11 (A; = 2). The 
post-break optical slope, and post-hidden break X-ray slope, 
could then be caused by a circumburst density profile of 
k — 1.3 ± 0.7. The pre-break slopes, which arc shallower 
than predicted by Aa = 0.97 ± 0.05, could be caused by 
continued energy injection with q = 0.8 ± 0.2, though this 
would increase the energy of the afterglow over the observed 
period of injection by a factor of / ~ 30. Given that we can 
equally well explain a cooling break either above or below 
the observed frequencies, we cannot say which of the two 
explanations is favoured. 



4.2.8 GRB 060729 

We find that p = 2.22 ±0.04 for < fx and q = 1.00 ±0.03 
with k unconstrained. Despite the fact that the broadband 
SED seems to favour a broken power-law (though this is 
not conclusive), the parallel nature of the lightcurves sug- 
gests that there is no cooling break between regimes so both 
the X-ray and optical regimes are either above, or below, 
the cooling break. For vx.opt > Vc, p = 2.22 ± 0.04 implies 
<a:x,opt,p = 1.17 ± 0.02 in agreement with both the X-ray 
and optical temporal slopes, ax, 2 and aopt,2. The break, to 
a = 1.70 ± 0.09, observed in the X-ray at lO" s is inter- 
preted as a jet break, albeit a slightly shallow one, though 
it may be rolling over to a steeper asymptotic value. Here, 
the early shallow phase is likely due to continued energy in- 
jection of the form, q = 1.00 ± 0.03, causing the observed 
X-ray break of Aa = 1.20 ± 0.05, and consistent with the 
optical. Such an energy injection would increase the energy 
of the afterglow over the observed period of injection by a 
factor of / ~ 10. 

If i^x.opt <Uc,p^ 3.22±0.04 implies ap = 1.67±0.02 or 
2.17±0.03 for a constant density or wind driven medium, re- 
spectively, steeper than either observed slope before ~ lO'' s 
but consistent with the slope after that time. The break with 
Aa > 0.39 ± 0.04, dependent on fc, can be explained by en- 
ergy injection of between q = 0.25±0.03 and q = 0.37±0.04. 
This would then mean that the early shallow phase can- 
not be explained by energy injection of the same form and 
would require an additional energy injection component, so 
we favour our initial explanation. 

4.2.9 GRB061121 

We find that p = 1.88 ± 0.08 for i^c < and A; = 1.5 ± 0.3, 
q = Q. The broken power-law SED implies that the X-ray 
regime is above the cooling frequency, which in turn implies 
p = 1.88 ± 0.08 and ax,p = 0.98 ± 0.04, in agreement with 
the observed X-ray index of ax,i = 0.97 ± 0.02 after 2000 s. 
The optical temporal slope, aopt.i, is then explained as being 
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Table 8. Most likely values for the blast wave parameters after our interpretation of individual bursts ('section l4.2| l. Where we interpret 
a jet break within the light curve, this is also noted. Values are bracketed where we cannot favour one interpretation over an other or 
the interpretation is not unequivocal. For those bursts with no valid blast wave interpretation, values are not included. 
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(-) 


(-) 
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I^X < 


2.64 
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-0.7 ± 1.4 


0.65 ± 0.15 
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i^X < 
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jet 
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(2.34 


± 0.16) 
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(jet) 
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(3.34 


± 0.16) 


(1.3 ± 0.7) 


(0.8 ± 0.2) 
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2.22 
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1.00 ± 0.03 


jet 
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Vc < i^x 


1.88 


± 0.08 


1.5 ± 0.3 





jet 


061126 















below the cooling break and in a medium with k = 1.8 ±0.3. 
The break observed in the X-ray light curve and consistent 
with the optical data is interpreted as a jet break rolling 
over to the expected value of ajct ~ p, in agreement with the 
observed value of Qx,2 = 1-55 ± 0.03. It is unclear what the 
X-ray emission before 2000 s is, but given that no shallowing 
is observed in optical, we can assume that it is not caused 
by continued energy injection, at least not in isolation. If 
it were due to continued energy injection, it would have a 
value of g ~ 0.55. 

4.2.10 GRB 061126 

We do not find a blast wave interpretation consistent with 
the light curves of this burst. The inconsistency of the opti- 
cal and X-ray slopes suggests that the cooling frequency is 
between the two regimes and p = 1.82 ± 0.08, in agree- 
ment with the SED which is best described with a bro- 
ken power-law. When i/x > Vc, the predicted value of 
Qx.p = 0.97 ± 0.04 significantly (8a") underestimates the ob- 
served value of ax = 1-29 ± 0.01, though the optical can be 
explained in this scenario if the circumburst density profile 
has an index of fe = 1.5 ± 0.3. The post-break optical slope 
of Q!opt = 1-9 ± 0.1 and an assumed 'hidden' X-ray break 
could be consistent with the predicted ajct = p in this inter- 
pretation. However, given the disagreement of the pre-break 
X-ray temporal slope with the predicted value, it seems that 
the data of this burst are incon sistent with th e blast wave 
model. This is consistent with iGomb oc et all (|2008l ) who 
suggest that the X-ray afterglow is due to an additional 
emission process rather than blast- wave shock emission. 



4.3 Distribution of p 

The universality of the electron energy dist ribution index, 
p, has been examined by several a uthors; Ch evalier & Li 
ll2000l).|Panaitescu fc Kumaij (|2002h , IShen et al.l (|2006 ) and 
IStarling et al.l 1 20081 ) applied different methods to samples 
of BeppoSAX and Swift bursts, all reaching the same con- 
clusion that the observed range of p values is not co nsis- 
tent with a single ce ntral value of p. IShen et al.l (|2006l ) and 
IStarling et al.l (|2008l ) showed that the width of the parent 
distribution is ~ 0.3 — 0.4. In the studies so far there have 
been some limitations: the BeppoSAX sample is limited. 



both in the number of GRBs and the temporal and spectral 
sampling; and the only stu dy of Swift bursts for this pur- 
pose so far by IShen et al.l only used the X-ray afterglows, 
which introduces a large uncertainty because the position of 
the cooling break is unknown. Here we examine the univer- 
satily of p for our sampl e of Swift bursts by u sing the same 
methods as described in IStarling" et all l|2008l ). We test the 
null hypothesis that the observed distribution of p (Table [H] 
excluding bracketed values and Figure ^ can be obtained 
from a parent distribution with a single central value of p. 

We minimise t he log-likelihood as given in equation 
3 of IStarling et al.l . and we find the most likely value of 
p = 2.21 ± 0.03. To test our null hypothesis, we generate 
10^ different synthetic data sets for p for the 6 bursts in 
our sample by taking random numbers from probability dis- 
tributions that are described by the most likely value of 
p = 2.21 and the 1 a uncertainties in the 6 values of p. We 
then take 6 values of p coming from these synthetic data 
sets and calculate the most likely value of p and the accom- 
panying log-likelihood. Comparing the log-likelihood of the 
synthetic data with that coming from the measured p values, 
we conclude th at the null hypo thesis is rejected at the > 5 a 
level. Following IStarling et all we can also derive the width 
of the parent distribution, which is 0.24 < o-gcat < 0.35 at 
the la level and 0.12 < o"scat < 1-22 at the 3o" level. We 
also tested the 10 values of p from Table [7] that offered the 
least deviation from the expected canonical value of ~ 2.2, 
with no regard for the interpretation of the light curves or 
SEDs. In this case, the most likely value of p is 2.19 ± 0.02 
and the values are inconsistent with one population at the 
5a level. 

This result confirms the results from previous studies 
and has important implications for theoretical particle ac- 
celeration studies. Some of these (semi-)analytical calcula- 
tions and simulations indi cate that there is a nearly universal 
value of p ~ 2.2-2.3 (e.g. lKirk et al.ll2000l : lAchterberg erall 
l200j ). while other studies suggest t hat there is a large range 
of possible values for p of 1.5 — 4 (|Baring||2004l ). Although 
we find that there is not a universal value of p, our val- 
ues for the width of the parent distribution indicate that 
it is not as wide as the latter study suggest. The width we 
find, (Tscat ~ .3 is comparab le to the numbers found by 
IShen et al.l and IStarling et all but our results are based on 
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Figure 2. Most likely values for p after our interpretation of 
individual bursts (Tabled; from left to right: GRB 050801, 
GRB 050802, GRB051109A, GRB 060124, GRB 060729, 
GRB 061121. The line represents the most likely value of p 
over the plotted sample (Section l4.3l l. 



a sample with better temporal and spectral sampling per 
GRB, on average. 

Given that the value of p is not universal, it may be pos- 
sible that it changes suddenly or evolves gradually with time 
or radius even in a single event as environmental blast wave 
parameters (e.g., m agnetic field, ambient density) change or 
evolve (e.g., Hamilton fc Petrosian 1992 : Vlahos ct al. 2004; 
iKaiseil 12003 ). If this is the case our derived values of p 
should be considered time averaged values of the parame- 
ter. A change or evolution of p should be observable as a 
change or evolution of the synchrotron spectral index, /3(p) , 
but should also cause a deviation from the power or broken 
power- law decay indices (q(p), Table|6]) of the hght curves, 
though neither of these effects is observable in our data. 



4.4 Circumburst density profile, k 

The density structure, or profile, of the circumburst medium 
is generally given as p, or n (number density), oc where 
fe = is a constant density, or ISM-like, medium and k = 2 
is a wind-like medium. The value of k has important impli- 
cations for the study of progenitor models as the currently 
favoured model for long GRBs, involving the collapse of a 
massive, Wolf-Rayet star, is expected to have an associated 
strong stellar wind affecting the circumburst environment. 
Support for a stellar wind has been observed via absorp- 
tion lines with relative velocities in UV spectra in a num- 



ber o f cases, most notably GRB 02100 4 (e.g., Schaefer et al 



l2003l: IStarliiig et al.l [2OO5I ) . However. iPanaitescu fc Kumaij 



(I2OO2') found, from broadband modeling of 10 bursts, that 
a constant density medium was favoured in half the cases 
and in only one cas e was a wind driven profile necessary. 
IStarling et al.l l|2008l ) found, from a sample of 10 bursts, that 
4 were likely wind driven while 2 were certainly not consis- 
tent with a wind. 

In our sample of 10 bursts, only 6 have X-ray or optical 
light curves below the cooling break, Vc, where they are de- 
pendent on the circumburst density profile and in the case of 
GRB 060206 this interpretation is only one of two that can- 
not be distinguished. Of our calculated values of k (Table[8] 



Figure 3. Calculated values of k after our interpretation of 
individual bursts (Tabled; from left to right: GRB 050801, 
GRB 050802, GRB051109A, GRB 060124, (GRB 060206), 
GRB 061121. The lines represent k = (constant density) and 
fc = 2 (wind like) media (Section l4.4l l. For bursts where we cannot 
favour one interpretation over an other or the interpretation is 
not unequivocal, we have bracketed the GRB name. 



and Figure[31l, 2 are consistent (2cr) with both fe = 2 and 
fc = (GRB 050801 and GRB 060206); 2 are best described 
with a wind-like medium and are inconsistent (5cr) with 
fc = (GRB 050802 and GRB 061121); and 2 are consistent 
(2.1(t) with a constant density medium but inconsistent with 
fc = 2 (GRB051109A and G RB 060124). Our results are 
henc e consistent with th ose of (|Panaitescu fc Ku mar 200^ 
and l|Starling et al.| [200^) insofar as the sample requires both 
constant density and wind driven media to explain the ob- 
served broadband emission, as the sample is not consistent 
with only one of the two environments. 

It seems clear that the circumburst environments are 
unlikely to be drawn from only one of the ISM or wind- 
like profiles and, in fact, the wind-like expectation is only 
valid assuming a constant mass-loss rate and wind veloc- 
ity form the progenitor star; precarious assumptions for a 
star in its final death throws and given that the circumburst 
density profile may be affected by a number of processes. 
Firstly, a reverse shock is expected to propagate backwards 
through the wind as it meets the interstellar medium and 
would cause a constant density, shocked wind. Also, the un- 
shocked ISM itself has an constant density profile but its 
distance from the burst makes it unlikely that we would ob- 
serve the blast wave traversing the ISM. However, recent 
observations of Wolf-Rayet binaries suggests that this shock 
front might be much closer to the progenitor if the wind is 
asymmetric l|Eldridgdl2007l ). This would make it more likely 
that the blast wave would already be travelling through the 
ISM at the time of observations, as would a high progenitor 
velocity towards the observer. In fact, the transition of the 
blast wave across the wind's shock front to the unshocked 
ISM has been suggest ed in a small number of cases, for in- 
stance, GRB 030329 l|Pe'er fc Wiierj |2006| ). Not only can 
we expect density profiles of fc = 2 or fc = 0; IStarling et al.l 
discuss the fact that in Type Ib/c supernovae, with which 
GRBs are closely associated (e.g., GRB 030329/SN 2003dh, 
iHiorth et al.ll200^ ). densi ty profiles are found to range from 
fc = 1.5tofc = 2.0 (see IWeiler et al.l [20041 and references 
therein). 
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Figure 4. Calculated values of q, including those with q = 
on line, after our interpretation of individual bursts (Table[8ll; 
from left to right: GRB 050801, GRB 050802, GRB 051109A, 
ORB 060124, (GRB 060206), GRB 060729, GRB 061121 (Section 
I4.5|l . For bursts where we cannot favour one interpretation over an 
other or the interpretation is not unequivocal, we have bracketed 
the GRB name. 



4.5 Rate of continued energy injection, q 

iNousek et all ()2006h proposed continued energy injection as 
a process to explain tlie siiallower tlian expected decay in- 
dices observed in many Swift bursts. We assume tliis takes 
the form of i5 oc wliere q is dependent on tlie siiallowing 
effect of tlie injection, Aa = — Qobsorvcd , and on the value 
of p. The continued energy injection is either due to i) a dis- 
tribution of the Lorentz factors of the shells ejected from the 
central engine which causes those shells with lower Lorentz 
factor to catch up with the blast wave at a later time, or 
ii) continued activity of the central engine itself. The for- 
mer has a limit of g < 3 — fc, wher e k is the dens ity profile 
of the circumburst medium, which iNousek et al.l suggested 
could be used as a diagnostic to differentiate between the 
two sources of continued energy injection. 

For 7 of the 10 light curves in our sample we are able to 
estimate the level of energy injection, or to say that it is con- 
sistent with zero (Tableland Figure|4]), though in the case 
of GRB 060206 there are two equally valid interpretations, 
one requiring energy injection while the other does not. Ex- 
cluding GRB 060206, energy injection is only required in 3 
out of 6 bursts and each of these has an injection index, q, of 
between 0.65 and 1.0, corresponding to increases in the en- 
ergy of the afterglows over the observed periods of injection 
of a factor, / of between 5 and 10. These values of q and 
/ are consistent with both scenarios of energy injection a nd 
with previously calculated values of INousek et al.l (|2006h . 

Clearly a much larger sample of multi- wavelength after- 
glows is required to be able to constrain the s ample prop- 
erties of q. These results compare to those of lEvans et al.l 
(|2008l ) who find that the majority of XRT light curves with 
observed breaks require energy injection. However, those au- 
thors study a significantly different, much larger sample of 
bursts while we consider a sample of only the most well 
sampled. 



5 CONCLUSION 

Throughout this paper we have applied the bl ast wave model 
ees fc Meszarosiri992l : lMeszaros et al ] |l998l ). assuming on- 
axis viewing, a standard jet structure and no evolution of the 
microphysical parameters, to a selection of GRBs with well 
sampled temporal and spectral data. We attempted to con- 
strain three parameters of interest and to test the validity of 
the blast wave model for light curves observed by Swift, the 
complexity of which has made the identification of breaks 
and the interpretation of the blast wave more difficult than 
in the pre-Swift era. We find that the majority of the after- 
glows are well described within the frame work of the blast 
wave model and that the parameters derived are consistent 
with those values found by previous authors. 

After an inspection of all the Swtft X-ray light curves 
and a comparison with the optical data available in the lit- 
erature, we identified a sample of 10 bursts that had enough 
data to well constrain the optical and X-ray temporal in- 
dices, and the X-ray spectral indices, of the afterglows. We 
analysed these data, fitting power-law decays to the light 
curves in an attempt to identify possible breaks; taking into 
account the possibility of hidden breaks and the effect of 
host galaxy contribution to the optical, where possible. 

In the case of well constrained X-ray spectra, we show 
that p can be estimated from the X-ray spectra alone, with 
the caveat that it will give two possible values of p; these 
must be differentiated by either broadband SEDs or a com- 
parison with the predictions of the blast wave model, or 
both. We compare the observed light curves and spectra 
of each burst in our sample to those values expected from 
the blast wave model and find that we can successfully in- 
terpret all bursts, except two, within this frame work; these 
interpretations require 'hidden' X-ray breaks in GRB 060206 
and GRB 061126. Furthermore, we identify, reasonably un- 
ambiguously, jet breaks in 5 afterglows out of our sample 
of 10, including a jet break in GRB 060729 that was not 
identified in previous analysis of the light curve. 

After interpretation within the blast wave model, we are 
able to confidently estimate the electron energy distribution 
index, p, for 6 of the bursts in our sample. A statistical anal- 
ysis of the distribution of p reveals that, even in the most 
conservative case of least scatter, the values are not consis- 
tent with a single, universal value suggested by some studies; 
this has important implications for theoretical particle ac- 
celeration studies. In a number of cases, we are also able to 
obtain values for the circumburst density profile index, k, 
and the index of continued energy injection, q. The calcu- 
lated values of q are consistent with both suggested sources 
of continued energy injection and a much larger sample of 
afterglows will be required to constrain the sample proper- 
ties of q. The values of k, consistent with previous works on 
the matter, suggest that the circumburst density profiles are 
not drawn from only one of the constant density or wind-like 
media populations. 
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